Absence of boron aggregates in superconducting silicon confirmed by atom 
probe tomography 

K. Hoummada, 1 F. Dahlem, 2 -^ T. Kociniewski, 3 J. Boulmer, 3 C. Dubois, 4 G. Prudon, 4 
E. Bustarret, 2 H. Courtois, 2 D. Debarre, 3 and D. Mangelinck 1 

IM2NP, CNRS and Universite Aix- Marseille, Marseille, France. 

2 ) Institut Neel, CNRS and Universite Joseph Fourier, Grenoble, France. 

3 ) Institut d'Electronique Fondamentale, CNRS and Universite Paris Sud, Orsay, France. 
^Institut des Nanotechnologies de Lyon, CNRS and INS A Lyon, Villeurbanne, France. 

Superconducting boron-doped silicon films prepared by gas immersion laser doping (GILD) technique are 
analyzed by atom probe tomography. The resulting three-dimensional chemical composition reveals that 
boron atoms are incorporated into crystalline silicon in the atomic percent concentration range, well above 
their solubility limit, without creating clusters or precipitates at the atomic scale. The boron spatial distri- 
bution is found to be compatible with local density of states measurements performed by scanning tunneling 
spectroscopy. These results, combined with the observations of very low impurity level and of a sharp two- 
dimensional interface between doped and undoped regions show, that the Si:B material obtained by GILD is 
a well-defined random substitutional alloy endowed with promising superconducting properties. 



In the last few years, gas immersion laser doping 
(GILD) technique has succeeded to provide ultra-shallow 
silicon junctions doped with boron atoms at a concentra- 
tion well into the percent range. Samples prepared by 
this out-of-equilibrium technique present metallic con- 
ducting properties at room temperature 1 and a super- 
conducting state below about 0.6 Kelvin 2 ' 3 . However, 
the incorporation of doping atoms well above the solubil- 
ity limit is usually thermodynamically unstable, leading 
to phase separation and formation of aggregates 4 . Such 
material modifications would drastically affect electrical 
properties, thus limiting potential applications. The re- 
cent observation of a spatial dispersion in the supercon- 
ducting energy gap of GILD silicon 5 raises specifically 
the important question of the boron spatial distribution 
in these epilayers. 

In this Letter, we report on atom probe tomography 
(APT) measurements of the boron distribution in heavily 
doped silicon epilayers prepared by GILD. This accurate 
three-dimensional spatial study shows that the boron 
atoms incorporated well above their solubility limit are 
still randomly distributed into the silicon lattice, without 
forming any cluster or precipitate. In such a well-defined 
random substitutional alloy, the residual inhomogeneity 
of the boron atoms density is discussed as the origin of 
the spatial evolution of the local density of states. 

All silicon samples in this study were doped by GILD 
under the same conditions with a few 10 21 boron atoms 
per cubic centimeter 1 . Starting from a 10 -9 mbar resid- 
ual pressure, a precursor gas atmosphere of boron trichlo- 
ride was injected into the growth vacuum chamber to sat- 
urate the chemisorption sites of a <100>-oriented pure 
silicon wafer surface. This surface was melted by a ho- 
mogenized 20 ns-long xenon chloride excimer laser pulse. 



Boron atoms then diffused into the silicon liquid, which 
re-crystallized from the Si/Si:B interface to the surface, 
in epitaxy with the substrate, see X-ray diffraction mea- 
surements in Ref. [2] or [3]. Such a boron injection/laser 
melting/re-crystallization cycle was repeated 500 times 
in order to heavily dope a 2 x 2 mm 2 area over an ex- 
pected thickness of 120 nm. The final doping level in the 
silicon epilayer yielded a related critical superconduct- 
ing temperature of 510 mK, see transport data set B in 
Ref. [5]. 

Atom probe tomography is a 3D compositional profil- 
ing analysis based on field evaporation of atoms 6 . In the 
present study, GILD Si/Si:B thin films were first shaped 
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FIG. 1. Reconstruction of a 60x60x200 nm 3 volume in heavily 
boron-doped silicon epilayers analyzed by atom probe tomog- 
raphy. Each dot corresponds to an atom. Representative of 
the whole distribution, only 1% (10%) of the captured silicon 
(boron) atoms are shown. A sharp interface is clearly visible 
between undoped and doped regions. 
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into tips with 50 nm typical diameter using the lift-out 
technique 7 and the annular milling method 8 . Pulsed 
picosecond-laser evaporation was then performed on in- 
dividual tips with an Imago LEAP 3000X HR APT in- 
strument equipped with a large-angle reflectron. The 
analyses were done at a specimen temperature of 50 K, 
a 100 kHz pulse repetition rate, an energy of 0.7 nJ (or 
0.5 nJ) per pulse, an evaporation rate of 0.01 ions per 
pulse and a pressure of less than 10 -11 mbar. In total, 
four tips made from two twin samples, i.e. samples pre- 
pared under the same GILD conditions, were analyzed. 
Similar results have been obtained for all tips. 

Fig. 1 shows a reconstruction of a 60 x 60 x 200 nm 3 
tip volume, which represents a typical result of our APT 
measurements. Because of the high efficiency of ana- 
lyzed atoms by APT (around 40% compared to few per- 
cent with secondary ion mass spectroscopy, SIMS), only a 
fraction of boron and silicon detected atoms is displayed 
in Fig. 1 as red and grey dots, respectively. In this mea- 
surement, pure and doped silicon regions are easily iden- 
tified from their composition. In the pure silicon region, 
no peak was found for boron atoms in the local mass spec- 
trum. The detection limit was ~ 5.10 18 at. /cm 3 . The 
interface between Si and Si:B is well described by an iso- 
concentration of 2.5% of boron atoms. It was found flat, 
as expected for a liquid phase epitaxy from a crystalline 
interface 1 , with a roughness of about 0.7 nm. Such a 2D 
sharp transition can be related to the much larger diffu- 
sion length of boron atoms in the liquid layer compared 
to the solid material. At the Si/Si:B junction, carbon, 
the dominant impurity, reaches a maximum concentra- 
tion at 0.25%, against about 0.05% throughout the film. 
No chlorine atoms brought by the dopant gas 9 was de- 
tected in the doped layer, confirming the purity of the 
film. 

Fig. 1 data provide information about Si:B region com- 
position at the nanometer scale, showing the important 
result that nanosized precipitates or clusters 4 of boron 
atoms are absent. In order to determine precisely the 
one-dimensional boron depth distribution, a radial in- 
tegration was performed over a 10 nm diameter cylin- 
drical region perpendicular to the Si/Si:B interface, see 
Fig. 2a. The related concentration profile as a function 
of the depth (square symbols) presents two distinct zones 
with an average boron atom concentration of about 10% 
(zone 1) and 8% (zone 2). These concentration levels are 
compatible with a boron incorporation rate estimated at 
1-1.2 10 14 at. /cm 2 per laser shot of 1 J/cm 2 [10]. Pre- 
vious works 2 ' 3 showed a clear difference of behavior in 
terms of SIMS profile, diffraction peaks and critical tem- 
perature between samples having concentration of boron 
atoms above or below a value of about 4 %. Occuring 
only in the case of very heavy doping, the presence of 
two zones result either from a crystalline relaxation or 
from the lower boron atomic mobility in a highly concen- 
trated boron silicon liquid. 

High concentration SIMS profiles measured in a twin 
sample (full line in Fig. 2a) have been used to calibrate 
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FIG. 2. (a) One-dimensional boron depth distribution (red 
square) calculated by integrating APT data of Fig. 1 over a 
10 nm radius cylindrical volume (see left drawing). The blue 
line corresponds to secondary ion mass spectroscopy measured 
with a CAMECA IMS4F on a twin sample. Vertical dotted 
line delimits the transient surface region, (b) Histogram show- 
ing the distribution of the boron contents in sets of 200 atoms 
located in 30 x 30 x 10 nm 3 data slice (position A) compared 
to a binomial shape (black line). 



the z-axis of the APT reconstruction and furthermore 
coincide with the APT curve. The weak discrepancy is 
attributed to the existence of slightly different boron con- 
centration between twin samples. Large deviations found 
near the surface (0 < z < 5 nm) are well-known measure- 
ment artefacts 11 . In principle, it is difficult to quantify 
high concentration levels above 1 at. % with SIMS due 
to matrix effects. The present quantitative similarity be- 
tween the APT and SIMS data confirms the accuracy of 
the recently proposed Isotopic Comparative Method used 
to calibrate SIMS measurements to the alloy concentra- 
tion range 12 . 

Using a statistical analysis of the large volume recon- 
struction shown in Fig. 1, we can evaluate the dispersion 
in boron concentration values and determine whether or 
not aggregates exist at the atomic scale. For these pur- 
poses, sets of fixed numbers of atoms have been extracted 
from a 30 x 30 x 10 nm 3 data slice. They have been af- 
terward counted and classified according to their boron 
concentration. For instance, Fig. 2b shows the obtained 
distribution at the depth A for sets containing 200 atoms 
(reduced x 2 ~ 1-1)- Several sets sizes varying from 300 
to 50 atoms for one set have also been tested. The same 
binomial lineshape characteristic of a random distribu- 
tion 13 was obtained for all set sizes. By considering the 
smallest atomic set size, i.e. 50 atoms, and that the boron 
concentration in a set varies between 2 to 11%, we con- 
clude to the absence of boron clusters on a scale as small 
as 17 boron atoms (this value takes into account the de- 
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FIG. 3. (a) Schematic showing the position of the data boxes 
used to determine the lateral boron distribution, (b) Ob- 
tained distribution in a x-y plane of zone 1 at z = 30 nm. 
The z integration over 10 nm was performed perpendicular to 
the Si/Si:B interface. The resolution is fixed by the integra- 
tion over a 1 nm 2 area. 



tection efficiency of 37%). Such a statement is made 
possible by the large measured volume, which provides 
from one dataset statistical information on the distribu- 
tion of solute elements down to the atomic scale. Sev- 
eral other tests such as first nearest neighbor method 7 , 
radial distribution functions (B-B), or cluster identifica- 
tion routines 14 ' 15 were applied to different regions: no 
atomic cluster could be detected, either. In low boron 
concentration samples, Si:B composition is known to be 
difficult to obtain by APT due to a high evaporation field 
of boron compared to silicon 16 . Here, in heavily boron- 
doped silicon, we might benefit from the random boron 
distribution. 

In comparison, for boron atoms implanted in silicon, 
precipitates of SiB3 stable phases form at 900° C. The 
formation of metastable boron clusters was even observed 
at temperature as low as 600° C. These clusters were 
linked to the high density of defects generated by the 
implantation process, such as Si self interstitials and 
their associated larger defects 4 . During GILD out-of- 
equilibrium doping process, the temperature is also high 
enough to form a SiB3 precipitate in a stable phase. It is 
therefore surprising that boron atoms at concentrations 
much higher than their solubility limit in crystalline sili- 
con (approximately 0.2 at.% at 1000° C) do not precipi- 
tate. It is likely the high speed of recrystallization (typ- 
ically 3-4 m/s) 10 and the lack of implantation defects in 
the GILD process that hinder the precipitation of SiB3 
and the related boron clustering. Kinetics overcome here 
thermodynamics. 

From our APT measurements, we have also extracted 
the 2D-mapping boron atoms lateral distribution by 
making an integration in the z-direction from a thin box. 
Several z-thickness between 2 to 10 nm have been tested. 
The two boxes in Fig. 3a represent the 10 nm case. As 
shown in Fig. 3b, boron concentration was found to vary 
over a typical length of 8 nm between 6.4 and 9.1% in 
zone 2, and between 7.9 and 12.1% in zone 1. Concen- 
trations are less spread compared to Fig. 2b since we 
averaged over a larger thickness. The deviation from 



the average value is higher for the area containing more 
dopants, which argues for an even better homogeneity 
in samples prepared with a lower boron atoms concen- 
tration. Although important modifications of the mate- 
rial composition like clusters and aggregates are absent 
in heavily-doped GILD Si epilayers, fluctuation in the 
distribution of uncorrelated boron atoms remains on a 
nanometer scale. 

If these boron atoms are isolated substitutional atoms 
that do not form diboron or B-H pairs, such an inho- 
mogeneous distribution of uncorrelated boron atoms is 
likely to be at the origin of hole density fluctuations and 
disorder and thus determine its relevant spatial scale. 
Above the superconducting transition temperature, this 
leads to a very dirty metal behavior 2,10 . In the super- 
conducting state, the nanometer-scale puddle-like inho- 
mogeneities should increase localization of the supercon- 
ducting Cooper pairs and thus modulate the spatial dis- 
tribution of the superconductor gap. 

Actually, superconductivity spatial fluctuations were 
recently observed on twin films by scanning tunneling 
spectroscopies 5 in the subKelvin temperature range see 
Fig. 4a. This detailed measurement of local electronic 
density of states showed that the superconductivity fol- 
low the Bardeen-Cooper-Schrieffer model and that its en- 
ergy gap A(r) varies by up to about ±15% on a lateral 
spatial scale of about 10 nm and more. Fig. 4b illus- 
trates this inhomogeneity by presenting the energy gap 
value measured at different spatial positions on the sur- 
face. Every tested position showed a distinct supercon- 
ductivity gap. 

Our experimental observations confirm that the behav- 
iors of the boron atoms concentration and of the super- 
conductor energy gap are compatible both in terms of 
relative magnitude variation (on the order of 15%) and 
characteristic length scale (of about 10 nm). Neverthe- 
less, the observed characteristic length scale appears as 
significantly smaller than the calculated superconducting 
coherence length of about 50 nm [3] . The usual approach 
combining BCS theory and Drude model, leading to the 
latter value of the coherence length, thus appears not 
fully relevant in a very dirty metal like GILD heavily 
boron-doped silicon epilayers with an elastic mean free 
path of the order of few nanometers. 

In conclusion, boron atoms incorporated inside a sili- 
con lattice by immersion laser doping technique do not 
form precipitates even when their concentration is well 
above the solubility limit. The doping distribution on a 
nanometric scale is considered responsible for the local 
fluctuations of the superconducting properties in Si:B. 
The random boron distribution combined with a good 
purity material and the observation of an abrupt Si/Si:B 
interface demonstrate the feasibility of building innova- 
tive devices based on such a GILD-doped silicon. Here, 
atom probe tomography appears as a perfect tool to anal- 
yse the three dimensional chemical composition of super- 
conducting semiconductor 17 19 . In the future, this tech- 
nique could be used for several other covalent cubic semi- 
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FIG. 4. (a) Normalized differential conductances (dl/dV) of 
tunnel contacts on superconducting silicon at two different 
spatial positions (T sarnp i e = 110 mK). (b) Distribution of the 
superconducting energy gap values obtained at different lat- 
eral distance from an arbitrary reference point on two twin 
samples. 



conductors showing superconductvity like diamond 20 , sil- 
icon carbide 21 or germanium 21,22 . 
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